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ABSTRACT 

High time resolution of ultra wideband (UWB) signals facil- 
itates very precise positioning capabilities based on time-of- 
arrival (TOA) measurements. Although the theoretical lower 
bound for TOA estimation can be achieved by the maximum 
likelihood principle, it is impractical due to the need for ex- 
tremely high-rate sampling and the presence of large number 
of multipath components. On the other hand, the conven- 
tional correlation-based algorithm, which serially searches 
possible signal delays, takes a very long time to estimate the 
TOA of a received UWB signal. Moreover, the first signal 
path does not always have the strongest correlation output. 
Therefore, first path detection algorithms need to be consid- 
ered. In this paper, a data-aided two-step TOA estimation al- 
gorithm is proposed. In order to speed up the estimation pro- 
cess, the first step estimates the rough TOA of the received 
signal based on received signal energy. Then, in the sec- 
ond step, the arrival time of the first signal path is estimated 
by considering a hypothesis testing approach. The proposed 
scheme uses low-rate correlation outputs, and is able to per- 
form accurate TOA estimation in reasonable time intervals. 
The simulation results are presented to analyze the perfor- 
mance of the estimator. 

Index Terms — Ultra-wideband (UWB), impulse radio 
(IR), time of arrival (TOA) estimation, statistical change de- 
tection, method of moments (MM). 



1. INTRODUCTION 

Since the US Federal Communications Commission (FCC) 
approved the limited use of UWB technology 1 1 1, communi- 
cations systems that employ UWB signals have drawn con- 
siderable attention. A UWB signal is defined to be one that 
possesses an absolute bandwidth larger than 500MHz or a 
relative bandwidth larger than 20%. It can coexist with in- 
cumbent systems in the same frequency range due to its large 
spreading factor and low power spectral density. UWB tech- 
nology holds great promise for a variety of applications such 
as short-range high-speed data transmission and precise lo- 
cation estimation. 
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Commonly, impulse radio (IR) systems, which transmit 
very short pulses with a low duty cycle, are employed to im- 
plement UWB systems (|2|-|6|). In an IR system, a train of 
pulses is sent and information is usually conveyed by the po- 
sition or the polarity of the pulses, which correspond to Pulse 
Position Modulation (PPM) and Binary Phase Shift Keying 
(BPSK), respectively. In order to prevent catastrophic col- 
lisions among different users and thus provide robustness 
against multiple-access interference, each information sym- 
bol is represented by a sequence of pulses; the positions of 
the pulses within that sequence are determined by a pseudo- 
random time-hopping (TH) sequence specific to each user 

El 

The high time resolution of UWB signals facilitates very 
precise TOA estimation, as suggested by the Cramer-Rao 
lower bound. However, in practical systems, the challenge 
is to perform this estimation in a reasonable time interval us- 
ing low-rate correlation outputs. 

The conventional correlation-based TOA estimation al- 
gorithms are both suboptimal and require exhaustive search 
among thousands of bins, which results in very slow TOA 
estimation J7J. In order to speed up the process, differ- 
ent search strategies, such as random search or bit reversal 
search, are proposed in 0. In 0, a generalized maximum 
likelihood (GML) estimation principle is employed to obtain 
iterative solutions after some simplifications. However, this 
approach requires very high rate sampling, which is not prac- 
tical in many applications. An alternative to the GML-based 
approach is a low complexity non-data-aided timing offset 
estimation technique based on symbol-rate samples based 
on the novel idea of "dirty templates" 1 1 1 - [ 13 1. The main 
disadvantage of the timing with dirty templates (TDT) algo- 
rithm is that its TOA estimate will have an ambiguity equal 
to the extent of the noise-only region between consecutive 
symbols. 

In this paper, we propose a two-step TOA estimation al- 
gorithm. In order to speed up the estimation process, the first 
step estimates the rough TOA of the received signal based on 
received signal energy. Then, in the second step, the arrival 
time of the first signal path is estimated from low-rate corre- 
lation outputs by considering a hypothesis testing approach. 

The remainder of the paper is organized as follows. Sec- 
tion II describes the transmitted and received signal models 
in a frequency-selective environment. The two-step TOA es- 
timation algorithm is considered in Section EH, which is fol- 
lowed by simulation results in Section IV, and concluding 
remarks are made in the last section. 



2. SIGNAL MODEL 



3. TWO-STEP TOA ESTIMATION ALGORITHM 



We consider a BPSK TH-IR transmitted signal represented 
by: 

OO 

stx(f) = ^f~E ^2 a J b [j/N f \Wtx(t - jT f - CjT c ), (1) 

3=— oo 

where iUt x (i) is the transmitted UWB pulse with duration 
T c , E is the transmitted pulse energy, Tf is the "frame" time, 
and b|j/jVfJ — 1} i s tne bi narv information symbol. 

In order to smooth the power spectrum of the transmitted 
signal and allow the channel to be shared by many users 
without causing catastrophic collisions, a time-hopping (TH) 
sequence cj G {0, 1, N c — 1} is assigned to each user, 
where N c is the number of chips per frame interval; that 
is, N c = Tf/T c . Additionally, random polarity codes, a/s, 
can be employed, which are binary random variables taking 
on the values ±1 with equal probability, and are known to 
the receiver. Use of random polarity codes helps reduce the 
spectral lines in the power spectral density of the transmitted 
signal 1 14 1 and mitigate the effects of MAI \ \T5i . 
Consider the following channel model 

L 

h{t) = aiS(t -(I- 1)T C - ttoa), (2) 
(=i 

where a; is the channel coefficient for the Ith path, L is the 
number of multipath components, and ttoa is the TOA of 
the incoming signal. 

From Q and (|2ji, we can express the received signal as 

L 

r{t) = V^a ; s rx (t -{I- 1)T C - ttoa) + n(i), (3) 
where s rx (£) is given by 

oo 

Srx(i) = a j b U/N f iW IX (t - jT f - c 3 T c ), (4) 

j=—oo 

with w IX (t) denoting the received UWB pulse with unit en- 
ergy. 

Assuming a data aided TOA estimation scheme, we con- 
sider a training sequence of bj = 1 Vj. Then, (0) can be 
expressed as 

oo 

Srx(i) = ^ a 3 w rx(t - jTf ~ CjT c ). (5) 
j=-oo 

We assume, for simplicity, that the signal always arrives 
in one frame duration (ttoa < Tf), and there is no inter- 
frame interference (IFI); that is, Tf > (L + c max )T c (equiv- 
alently, N c > L + c max ), where c max is the maximum value 
of the TH sequence. Note that the assumption of ttoa < Tf 
does not restrict the validity of the algorithm. In fact, it is 
enough to have ttoa < T s for the algorithm to work when 
the frame is large enough and predetermined TH codes are 
employed. Moreover, even if ttoa > T s , an initial energy 
detection can be used to determine the arrival time within a 
symbol uncertainty before running the proposed algorithm. 



One of the most challenging issues in UWB TOA estima- 
tion is to obtain a reliable estimate in a reasonable time inter- 
val under the constraint of sampling rate. In order to have a 
low power and low complexity receiver, one should assume 
low sampling rates at the output of the correlators. How- 
ever, when low rate samples are employed, the TOA esti- 
mation can take a very long time. Therefore, we propose 
a novel two-step TOA estimation algorithm that can perform 
TOA estimation from low rate samples (typically on the order 
of hundreds times slower sampling rate than chip-rate sam- 
pling) in a reasonable time interval. In order to speed up the 
estimation process, the first step estimates the coarse TOA of 
the received signal based on received signal energy. Then, 
in the second step, the arrival time of the first signal path is 
estimated by considering a hypothesis testing approach. 
Express the TOA ttoa in as 3 

ttoa = kT c = k b T b + k c T Cl (6) 

where k <E [0, N c — 1] is the TOA in terms of the chip interval 
T c , is the block interval consisting of B chips (TJ, = BT C ), 
and fcb G [0, N c /B — 1] and k c 6 [0, B — 1] are the integers 
that determine, respectively, in which block and chip the first 
signal path arrives. 

The two-step TOA algorithm first estimates the block in 
which the first signal path exists; then it estimates the chip 
position in which the first path resides. In other words, it can 
be summarized as: 

• Estimate kb from received signal strength (RSS) mea- 
surements. 

• Estimate k c (equivalently, k) from low rate correlation 
outputs using a hypothesis testing approach. 

3.1. First Step: Coarse TOA Estimation from RSS Mea- 
surements 

In the first step, the aim is to detect the coarse arrival time 
of the signal in the frame interval. Assume, without loss of 
generality, that the frame time Tf is an integer multiple of 
Tf,, the block size of the algorithm; that is, Tf = N b T b . 

In order to have reliable decision variables in this step, 
energy is combined from N\ different frames of the incom- 
ing signal for each block. Hence, the decision variables are 
expressed as 

Ni-l 

Y i=J2 F< .i' (7) 

3=0 

for i = 0, . . . , N b — 1, where 

pjT f +(i+l)T b +CjT c 

Y itj = / \r(t)\ 2 dt. (8) 

JjTf+in+CjTc 

Then, k b in is estimated as 

k b = arg max Yt. (9) 

0<i<N b -l 

3 For simplicity, the TOA is assumed to be an integer multiple of the chip 
duration T c . In a practical scenario, sub-chip resolution can be obtained 
by employing a delay-lock-loop (DLL) after the TOA estimation with chip- 
level uncertainty 1161 . 



In other words, the block with the largest signal energy is 
selected. 

The parameters of this step that should be selected appro- 
priately are the block size Tf, (Nb) and the number of frames 
N%, from which energy is collected. 

3.2. Second Step: Fine TOA Estimation from Low Rate 
Correlation Outputs 

After determining the coarse arrival time from the first step, 
the second step tries to estimate k c in Ideally, k c £ 
[0, B — 1] needs to be searched for TOA estimation, which 

corresponds to searching k £ [kbB, (kb + \)B — 1], with kb 
obtained from However, in some cases, the first signal 
path can reside in one of the blocks prior to the strongest one 
due to multipath effects. Therefore, instead of searching a 
single block, k £ [k b B - M u (k b + l)B - 1], with Mi > 0, 
can be searched for the TOA in order to increase the proba- 
bility of detection of the first path. In other words, in addition 
to the block with the largest signal energy, we can perform an 
additional backwards search over M\ chips. For notational 
simplicity, let U = {n s , n s + 1, . . . , n e } denote the uncer- 
tainty region, where n s = k\B — Mi andrt e = (kt, + l)B — 1 
are the start and end points. 

In order to estimate the TOA with chip-level resolution, 
we consider correlations of the received signal with shifted 
versions of a template signal. For delay iT c , we obtain the 
following correlation output 



iT c +N 2 T f 



r(t)Stemp (t - iT c ) dt, 



(10) 



IT,. 



where iV2 is the number of frames over which the correlation 
output is obtained, and s temp (t) is the template signal given 
by 

ip(*) = X! a 3 w ^( t ~ 3 T f - c i T c)- (11) 

3=0 



Sic 



From the correlation outputs for different delays, the aim 
is to determine the chip in which the first signal path has ar- 
rived. By appropriate choice of the block interval T b and Mi, 
and considering a large number of multipath components in 
the received signal, which is typical for indoor UWB sys- 
tems, we can assume that the block starts with a number of 
chips with noise-only components and the remaining ones 
with signal plus noise components, as shown in Figure [2 
Assuming that the statistics of the signal paths do not change 
significantly in the uncertainty region U, we can express the 
different hypotheses approximately as follows: 



H 
7~Lk 



for k £ U, where rji's denote the i.i.d. output noise dis- 
tributed as j\f(0 , cT^) with a 2 n = N 2 jV /2, a\, . . . , a nf -k+i 
are independent channel coefficients, assuming nf—n s + l < 
L, and rif = n e + with A/2 being the number of correla- 
tion outputs that are considered out of the uncertainty region 
in order to have reliable estimates of the unknown parameters 
of a. 
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Fig. 1. Illustration of the two-step TOA estimation algo- 
rithm. The signal on the top is the received signal in one 
frame. The first step checks the signal energy in Nb blocks 
and chooses the one with the highest energy (Although one 
frame is shown in the figure, energy from different frames 
can be collected for reliable decisions). Assuming that the 
third block has the highest energy, the second step focuses 
on this block (or an extension of that) to estimate the TOA. 
The zoomed version of the signal in the third block is shown 
on the bottom. 



Due to very high resolution of UWB signals, it is appro- 
priate to model the channel coefficients approximately as 



Oil 



di|ai|, 
di\ai\ 
0, 



P 
1 



1 = 2, 



(13) 



where p is the probability that a channel tap arrives in a given 
chip, di is the sign of a/, which is ±1 with equal probabil- 
ity, and I a; I is the amplitude of ai, which is modelled as a 
Nakagami-771 distributed random variable with parameter ft; 
that is fTTl . 



p(a) 



T(m) 



my 

n) 



„2m-l. 



(14) 



for a > 0, m > 0.5 and il > 0, where T is the Gamma 
function 1181 . 

From the formulation in H2\ . it is observed that the TOA 
estimation problem can be considered as a "change detec- 
tion" problem 1191 . Let 6 denote the unknown parameters 
of the distribution of a; that is, 6 = \p m O]. Then, the 



log-likelihood ratio (LLR) is calculated as 

Pe{zi\Hk) 
p{zi\Ho) ' 



ft f 

s n k f(8) = yiog 



i=k 



(15) 



where pe{zi\Hk) denotes the probability density function 
(p.d.f.) of the correlation output under hypothesis Hk and 
with unknown parameters given by 9, and p(zi\Ho) denotes 
the p.d.f. of the correlation output under hypothesis Ho- 

Since 9 is unknown, its ML estimate can be obtained first 
for a given hypothesis Hk and then that estimate can be used 
in the LLR expression. In other words, the generalized LLR 
approach 1 19 1 can be taken, where the TOA estimate is ex- 
pressed as 



where 



k = argmax S^ f (9 ML (k)), 

k£U 



0ML(fc) =arg sup 
e 



(16) 



(17) 



However, the ML estimate is usually very complex to cal- 
culate. Therefore, simpler estimators such as the method of 
moments (MM) estimator can be employed to obtain those 
parameters. The nth moment of a random variable X having 
Nakagami-m distribution with parameter fl is given by 



E{X n } = 



T{m + n/2) ( 0, 



T(m) 



rn 



n/2 



(18) 



Then, from the correlator outputs {zi}"i fc+1 , the MM es- 
timates for the unknown parameters can be obtained after 
some manipulation: 



7i 72 
272 - 73 



, rajfM 



272 ~ 73 



ft i 



9 , » l MM 
73 - 72 



27 2 2 ~ 73 

72 

(19) 



where 



Let pi(z) and P2{z), respectively, denote the distribu- 
tions of r\ and N2VE d\a\ + ?y. Then, the generalized LLR 
for the fcth hypothesis is given by 



,P2{Zk) 



n f 



tog £ 1 0| 



PP2{zi) + (l-p)pi(zi) 



' Pi(z k ) 



i=k+l 



where 



and 



with 



Pi(z) 



P2(z) 



V! —4- / 1 Z 1 

= — e 2 "™<P to, -; — 

2-Ka n V 2 V2 



(23) 



(24) 



(25) 



A 2V^r(2m) 

~ r(TO)r(TO + o.5) 



u 2 = 2cr„ 1 + 2to 



(26) 



and $ denoting a confluent hypergeometric function given 
by GU 



ftl! /3 2 (/3 2 + l)2! 
/3 1 (/3 1 + l)(/3 1 + 2) x 3 
/3 2 (/? 2 + l){/3 2 + 2) 3! 



(27) 



Note that the p.d.f. of N 2 v E d\a\ +rj, P2{z), is obtained 
from dl4> . (I24> and the fact that rfis ±1 with equal probabil- 
ity. 

After some manipulation, the TOA estimation rule can be 
expressed as 
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7i 
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2„4 
n / ? 

(20) 



with //j denoting the jth sample moment given by 

n f 
1 <5 

Mi 



^ zl 



nt — k ^— ' 

i=fe+l 



(21) 



Then, the index of the chip having the first signal path 
can be obtained as 



k = argmax S^ f (9 M M(k)), 
keu 



(22) 



where 8 M M(k) = [pmm m m m !1mm] is the MM estimate 
for the unknown parameters 4 . 

4 Note that the dependence of Pmm> m MM ar| d ^mm on the change 
position k is not shown explicitly for notational simplicity. 



k = arg max< loe 

keu 
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p to, 0.5 



V-2 



1-p 



(28) 



3.3. Additional Tests 

Note that the formulation in d!2t assumes that the block al- 
ways starts with noise-only components, and then the sig- 
nal paths start to arrive. However, in practice, there can be 
cases where the first step chooses a block consisting of all 
noise components. By combining a large number of frames; 
that is, by choosing a large Ni in 0, the probability of this 
event can be reduced considerably. However, very large Ni 
also increases the estimation time; hence, there is a trade- 
off between the estimation error and the estimation time. In 
order to prevent erroneous TOA estimation when a noise- 
only block is chosen, a one-sided test can be applied us- 
ing the known distribution of the noise outputs. Since the 
noise outputs have a Gaussian distribution, the test reduces 
to comparing the average energy of the outputs after the es- 
timated change instant to a threshold. In other words, if 



\ — y nf ~ zf < 61, the block is considered as a noise- 

nf — k+l £ -~ i %=K 1 

only block and the two-step algorithm is run again. 

Another improvement of the algorithm can be obtained 
by checking if the block consists of all signal paths; that 
is, the TOA is prior to the current block. Again, by fol- 
lowing a one-sided test approach, we can check the aver- 
age energy of the correlation outputs before the estimated 
TOA against a threshold and detect an all-signal block if the 
threshold is exceeded. However, for very small values of the 
TOA estimate k, there can be a significant probability that 
the first signal path arrives before the current observation re- 
gion since the distribution of the correlation output after the 
first path includes both the noise distribution and the signal 
plus noise distribution with some probabilities as shown in 
(II 31 . Hence, the test may not work although the block is an 
all-signal block. Therefore, some additional correlation out- 
puts before k can be employed as well when calculating the 
average power before the TOA estimate. In other words, if 

7 1 tz y\jZl t f zf > 82, the block is considered as an 

all-signal block, where M3 > additional outputs are used 
depending on k. When it is determined that the block con- 
sists of all signal outputs, the TOA is expected to be in one 
of the previous blocks. Therefore, the uncertainty region is 
shifted backwards, and the change detection algorithm is re- 
peated. 

4. SIMULATION RESULTS 




Fig. 2. Root mean square error (RMSE) versus signal-to- 
noise ratio (SNR) for four different IEEE 802. 15.4a channel 
models. The averages over 100 realizations are plotted for 
each model. 



tions. More accurate results can be obtained by employing 
different parameters in different scenarios. 



In this section, we perform simulations in order to eval- 
uate the performance of the estimator over realistic IEEE 
802.15.4a channel models 1 17 1. We consider residential and 
office environments with both line-of-sight (LOS) and non- 
line-of-sight (NLOS) situations. 

In our simulation scenario, the signal bandwidth is 
7.5GHz and the frame time of the transmitted training se- 
quence 300ns. Hence, we consider an uncertainty region of 
2250 chips, and divide the region into Nj, — 50 blocks. The 
number of pulses, over which the correlations are taken in 
the first and second steps is given by Ni — 50 and N2 — 25, 
respectively. Also Mi = 180 additional chips prior to the 
uncertainty region determined by the first step are included 
in the second step. 

The estimator is assumed to have 10 parallel correlators 
for the second step. In a practical setting, the estimator can 
use the correlators of a RAKE receiver that is already present 
for the signal demodulation, and 10 is a conservative value in 
this sense. 

In Figure|2] root mean square errors (RMSEs) are plotted 
for different signal-to-noise ratios (SNRs) for four different 
channel models. When the SNR is high enough, very ac- 
curate TOA estimates can be obtained. Due to the different 
characteristics of the channels in residential and office envi- 
ronments, the estimates are better in the office environment. 
Namely, the delay spread is smaller in the channel models for 
the office environment. Moreover, as expected, the NLOS 
situations cause increase in the RMSE values. 

Each TOA estimation is performed in approximately 0.92 
millisecond. Because we did not employ any additional tests 
afte r the TOA estimate, which are described in Sub-section 
13.31 and used the same parameters for all the channel models, 
the estimation time is the same for all the channel realiza- 



5. CONCLUSIONS 

In this paper, we have proposed a two-step TOA estimation 
algorithm, where the first step uses RSS measurements to 
quickly obtain a rough TOA estimate, and the second step 
uses a change detection approach to estimate the fine TOA 
of the signal. The proposed scheme relays on low-rate corre- 
lation outputs, but still obtains a considerably accurate TOA 
estimate in a reasonable time interval, which makes it very 
practical for UWB systems. 

The future work includes optimization of the estimator 
parameters, such as 2V&, by considering the statistics of the 
channel parameters defined by the IEEE 802.15.4a channel 
committee. 
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